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The formation of extended electron states in one-dimensional nanostructures is of key 
importance for the function of molecular electronics devices. Here we study the effects 
of strong electron—phonon interaction on the formation of extended electronic states in 
intentionally created Cl vacancy pairs and chains in a NaCl bilayer on Cu(lll). The in¬ 
teraction between the vacancies was tailored by fabricating vacancy pairs and chains of 
different orientation and separation with atomic precision using vertical manipulation. 

Small separation of divacancies led to the formation of symmetric and antisymmet¬ 
ric vacancy states and localized interface-states. By scanning tunneling spectroscopy 
(STS) we measured their energy splitting and broadening as a function of the inter¬ 
vacancy separation. Unexpectedly, the energy splitting between the vacancy states 
is enlarged by level repulsion resulting from phonon dressing of the electronic states, 
as evidenced by theory. Already for a few coupled vacancies we observe an emerging 
band structure of the defect band. 

One-dimensional nanostructures assembled from 
single metal atoms and molecules on surfaces are of 
great fundamental and technological interest as their 
electronic states exhibit a wealth of new quantum 
phenomena that might be exploited in future molecular 
electronics devicei^. Some examples of such phenomena 
identihed by scanning tunneling microscopy (STM) 
and spectroscopy (STS) include metal-atom chains 
exhibiting ’particle-in-box’ stateJ^, bistable antiferro- 
magnetisirP, Tomonaga-Luttinger liquid behavioui^^, 
and signs of Majorana spinorJ^. The formation of 
extended or band-like electronic states in metal-atom 
and molecular chains on surfaces has been widely studied 
experimentall}^^^^, but potentially interesting effects 
of strong electron-phonon coupling on these states, 
such as, for instance, polaronic effects, have attracted 
much less attention. A notable exception is the ob¬ 
servation of coherent electron-phonon states and the 
breakdown of the Born-Oppenheimer approximation in 
oligothiophene molecular chain^^. In the electronically 
adiabatic regime, where the Franck-Gondon principle 
is applicable, the electron-phonon coupling results in 
vibronic satellite structures. How these vibronic effects 
influence the delocalization of an injected electron in 
a one-dimensional nanostructure is an open question, 
in particular when the electron-phonon coupling is so 
strong that these effects dominate. 

Here we address this question in a combined ex¬ 
perimental and theoretical study of coupled electronic 
states localized at intentionally created G1 divacancies 


and vacancy chains in a NaGl bilayer on Gu(lll). The 
G1 vacancies provide a quantum-well structure that 
co uld b e created with atomic precision using the STM 
tip| i2 | i3 | energies of the electronic states formed 

in these structures by injection of an electron from 
the tip of an STM and their spatial behaviour on the 
atomic scale were resolved by STS and analyzed using 
simple models. Our findings have ramifications for the 
formation of extended states from coupled adsorbate 
or confined defect states (such as in dopanti^ or 
dangling bond^^ in semiconductors and ionic crystals. 
Hence, it could improve the understanding of polaronic 
effects in coupled quantum system s, such as quantum 
dot jlMUl, charge-based qubitJ^^^ and quantum cellular 
automateP^. 

Goupled, individual G1 vacancies in a NaGl bilayer 
on a Gu(lll) surface are ideal model systems to study 
the effect of strong electron-phonon coupling on the 
delocalization of electronic states. First, G1 monovacan¬ 
cies are atomically well-defined and highly symmetric 
defects. They feature two distinct electronic states: 
an unoccupied vacancy state (VS) that is strongly 
coupled to optical phonons in the film and a localized 
interfac e-stat e (LIS) with a negligible electron-phonon 
couplinj^^^. The electronic coupling between VSs or 
LISs of different vacancies can be tuned in a controlled 
way by the lateral spacing between vacancies with 
atomic precisiorP^. Finally, G1 vacancies feature a 
much higher stability against inelastic excitations than 
adsorbateJ^. G1 vacancies could even be filled with 
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other atomJ^. Hence, they allow us to explore artificial 
coupled quantum systems with great control. 

In a previous studj^, single Cl vacancies in a NaCl 
bilayer on copper surfaces were characterized in detail. 
A Cl vacancy in the outermost layer of this supported 
NaCl bilayer is analogous to the widely studied and 
well-known colour centers in bulk NaCP^ with one key 
difference; the localized VS is unoccupied because the 
electron can tunnel into unoccupied metal states. The 
unoccupied VS then gives rise to a positively charged 
vacancy and an attractive potential that is able to 
split off a LIS from the free-electron-like interface-state 
band of the NaCl/Cu(lllThis LIS was observed 
as a narrow resonance just below the bottom of the 
interface-state band in the differential conductance 
(d//dV) spectra, whereas the VS showed up as a broad 
negative ion resonance at a sample voltage of about 2.8 V 
in dl /dV spectra. The negative ion resonance was found 
to have a large Gaussian broadening (full-width half 
maximum (FWHM) of cr = 0.27V) because of the strong 
electron-phonon coupling of the electron in the VS to 
optical phonons in the NaCl bilayer. Despite this strong 
coupling corresponding to an estimated Huang-Rhys 
parameter of about S' = 30, the associated relaxation 
energy of the electron in the VS is not large enough 
to prohibit tunneling into unoccupied metal states and 
to allow the formation of a stable occupied VS state 
corresponding to a localized polaron. In contrast to 
the VS, the LIS was found to be dominated by lifetime 
broadening with a negligible phonon broadening in the 
d//dV spectra, owing to the efficient screening of the 
electron-phonon interaction by the metal electrons. 

The atomically precise locations of intentionally 
created Cl divacancies in a NaCl bilayer on Cu(lll) 
were resolved unambiguously using noncontact atomic 
forc e mic roscopy (AFM) with CO-terminated tips (CO 
tipsj33|34| Divacancies in the nineth-nearest-neighbour 
(9NN), 6NN, 5NN, 4NN, 3NN and 2NN configurations 
(see schematic in Fig. Ef) were created from two Cl 
vacancies close to each other but not in the vicinity 
of any other vacancies or defects. Constant-height 
AFM images of the 5NN to 3NN divacancies and of a 
single vacancy are shown in Fig. -e. The electronic 
structure of the VSs and the LISs of the divacancies were 
characterized by constant-height dI/dV{x^V) maps 
along the line connecting the two vacancies in Fig. [Ife-j 
and Fig. [^-o, respectively. All d//dV data shown are 
recorded using Cu-terminated tips. For the 4NN, 3NN 
and 2Nr'J^ divacancy configurations, two peaks with 
Gaussian line shapes were observed around 2.8 V. The 
dI/dV{x^V) maps show that the state at lower energy 
is localized between the vacancies (symmetric), whereas 
the state at higher energy has a nodal plane between 
them (antisymmetric), in analogy with the bonding 


and antibonding orbitals of the hydrogen molecule. 
Similarly, the divacancy LISs were also observed to 
split into a symmetric and an antisymmetric state. 
However, for small vacancy separations, the latter state 
lies above the NaCl/Cu(lll) interface band onset at 
Eb = —230 mV and overlaps with the band continuum. 
Hence, only the symmetric state is observed as a 
well-defined peak in d//dV. The splitting of both the 
VSs and LISs was found to increase with decreasing 
inter-vacancy distance. Tab. SI summarizes the observed 
peak positions and energy splittings of the VSs and LISs. 

To study the delocalization of VSs and LISs further, 
short one-dimensional arrays of N vacancies, referred to 
in the following as vacancy chains, were created in the 
5NN (apolar NaCl direction) and the 3NN (polar NaCl 
direction) configuration. They were also characterized 
by constant-height CO tip AFM images (Fig. i>.c). 
For the 5NN chains, the LISs form one-dimensional 
quantum-well states, observed as distinct resonances 
in d//dV, with an increasing integer number of nodal 
planes with voltage (see Fig. §> and Fig. S2). Sur¬ 
prisingly, also broad resonances are observed up to 
200 — 300 mV above the interface-state band onset. 
These quasi-bound LISs have a lower height in d//dV, 
owing to the increased lifetime broadening compared 
to the bound LISs. This originates from the additional 
decay channel into the interface-state band. Further¬ 
more, the energy position of the ground-state LIS does 
not change significantly from the V = 4 to the N = 6 
chain and seems to converge to about —0.33 V. The VSs 
in the 5NN chains only show a single resonance peak at 
about 2.7 V because of the weak electronic interaction 
between the VSs at these large intervacancy distances 
(cf. Fig. [^). In contrast, the 3NN chains displayed in 
Fig-i: have a significant interaction between the VSs 
and exhibit several resonances, as shown in Fig. The 
increasing number of nodal planes with voltage shows 
that the VSs are delocalized over these chains despite 
the strong electron-phonon coupling. In addition, the 
band width of VSs increases with chain length N, 
and is more or less saturated at 0.54 V for the V = 5 
chain. Interestingly, the level broadening decreases 
for longer chains as displayed in Fig. This finding 
can be rationalized by the larger spatial extent of the 
eigenstates, which decreases their relaxation energy and 
accordingly the effective electron-phonon coupling. 

In the following, the VS energies of the 3NN chains are 
compared with the results from a simple tight-binding 
(TB) model. First, we assume that each vacancy in 
the chain has an on-site energy eq and is coupled 
to its next neighbour(s) with a hopping term —t, as 
sketched in Fig. By adjusting the free fitting 

parameters Eq and t to the experimental mono- and 
divacancy VS energies, one can calculate the energy 
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levels for different chain lengths as shown in Fig. H- 
In comparison with the experimental energy spectrum, 
depicted in Fig. ii. the level spacings for each chain 
could be roughly reproduced, but the TB model does 
not capture the shift of the band center towards lower 
energies with increasing chain length. This shift can 
be explained by the electrostatic interaction between 
adjacent vacancies that deepens the potential at each 
vacancy and therefore lowers the on-site energy. To 
account for the energy shift we introduce the parameter 
S, which is the on-site energy shift due to a neigh¬ 
bouring vacancy (see Fig. §). The resulting energy 
levels of this extended TB model are shown in Fig. If. 
The calculated energy spectrum now agrees very well 
with the experimentally observed VS levels of the chains. 

So far, the strong electron-phonon coupling, as shown 
by the large broadening of the VSs, was neglected in 
our tight-binding models. The role of this coupling 
in the formation of delocalized states by a tunneling 
electron in vacancy pairs (and chains) needs to be 
investigated. Here, we use a simple generalization of 
the electron-phonon interaction model of the single 
vacanc}^ to vacancy pairs and chains. As detailed in 
the SI, the direct electronic interactions between the 
VSs are described by a simple tight-binding model, and 
each VS is coupled linearly to the phonon modes of the 
NaCl him. As we will show next, in this model, the 
strong electron-phonon coupling has a profound effect 
on the energy splitting and broadening of the vacancy 
states. This effect is illustrated here for a divacancy 
using appropriate parameters, where each single vacancy 
state is coupled to a single Einstein phonon mode. Here, 
the electronic interaction between the vacancies is large 
compared to the phonon energy, corresponding to an 
electronically adiabatic regime. 

In Fig. the adiabatic potential energy surfaces 
(PES) of the singly occupied, divacancy states for 
symmetric displacements Qs (where qi = ^ 2 ) of the two 
Einstein phonon modes are shown. They are split by 
2 t, where —t is the direct electronic interaction between 
the VSs, and both have the same shape as for the 
adiabatic PES of a monovacancy. In contrast, along 
antisymmetric displacements qas (where qi = — ^ 2 ) 
(Fig. 1^), the corresponding PESs exhibit an avoided 
crossing with a minimum energy gap of 2t. This level 
repulsion will be shown to have a profound influence on 
the phonon broadening and the formation of delocalized 
states. 

As for the single vacancy, the result for the phonon 
broadening of the tunneling through the divacancy 
states turns out to have a simple physical form in the 
prevailing strong electron-phonon coupling limit, as 
detailed in the SI. The phonon broadening is determined 


by probabilities for vertical Eranck-Condon transitions 
from the vibrational ground-state of the PES of the 
electronic ground-state to vibrational states on the 
two PESs for the divacancy states (gray arrows in 
Fig. If): and d). In the strong electron-phonon coupling 
limit, the transition probabilities are dominated by the 
contribution from the linear electron-phonon coupling 
terms. The corresponding phonon broadenings are then 
simply obtained from the change in the divacancy energy 
levels Eb (bonding state) and £„ (antibonding state) with 
the Gaussian fluctuations of the energies fei and Ss 2 of 
the vacancy 1 and 2 from their mean value Eq due to the 
zero-point motions of the phonons. As shown in the SI, 
this result is also valid for more general phonon baths 
and electron-phonon couplings. 

The strong coupling result for the phonon broadening 
of the divacancy states is illustrated in Eig. and 
d. A symmetric fluctuation Ses (where Sei = 5e2) 
results, as shown in Eig. |^, in a linear dependence of 
Eh and Ea on 6es with slope one, but they are split by 
2t. These Gaussian fluctuations then give rise to two 
overlapping Gaussian lineshapes centered at Eq ±t with 
a variance of {Se^) for the local density of states (LDOS). 

In contrast, the level repulsion results in a mini¬ 
mum energy separation 2t between Eh and Ea for an 
antisymmetric fluctuation (where Sei = —SE 2 ). 

Eurthermore, ^Eas breaks the symmetric and antisym¬ 
metric character of the two coupled VSs, and they 
tend to localize on each vacancy for SEas ^ t. The 
resulting LDOS from these Gaussian fluctuations, shown 
in Eig. B. consists of two relatively sharp, highly 
asymmetric lineshapes, with mean energies separated 
by more than the energy splitting 2t. In general, the 
fluctuations 5ei and 5e2 of the two vacancy energies 
will be a superposition of both 5es and SEas- In the 
case of uncorrelated energy fluctuations {5ei5e2) = 0 , 
the LDOS is a convolution of the lineshape from 5Eas 
with a Gaussian lineshape with a variance of {5e\ 2 ) /2 
because 5es and 5Eas contribute equally to (<^£ 12)5 
the variance of the single vacancy energy fluctuations. 
As the Gaussian shape dominates the broadening in 
such a convolution, the resulting LDOS as displayed 
in Eig. 1^ has a similar shape as two Gaussians with 
a EWHM reduced by almost l/\/2 as compared to an 
isolated vacancy. This decrease of the broadening in 
the uncorrelated case is the underlying reason why the 
observed broadenings of the levels of the chains decrease 
with chain size (see SI). 

The two phonon modes considered here affect only 
one of the VSs each and result in uncorrelated vacancy 
site energy fluctuations. In reality, however, a given 
individual phonon mode may act on the energy of 
both vacancies in a symmetric or antisymmetric fashion 
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(consider e.g. the motion of an ion centered between 
two vacancies). Hence, partially correlated vacancy site 
energy fluctuations cannot be excluded. As shown in 
the SI, the two Einstein modes and their symmetric and 
antisymmetric fluctuations can serve as a basis for the 
description of any phonon bath and electron-phonon 
coupling. Our experimental results are consistent with 
the assumption of fully uncorrelated fluctuations of the 
two VSs. Interestingly, in the uncorrelated case, the 
LDOS, shown in Fig. (red line), exhibits two peaks 
with an apparent separation that is substantially larger 
than the intrinsic splitting of 2t. This result is a char¬ 
acteristic effect of the level repulsion between ih and 
and is in sharp contrast compared to a superposition of 
two monovacancy resonances separated by the intrinsic 
splitting (shown in Fig. (blue line)). In the latter 
case, one would not be able to identify separate peaks 
due to the larger broadening and smaller splitting. 

The enhancement of the apparent peak separation 
in the uncorrelated case over the intrinsic splitting is 
shown in Fig. as a function of the scaled intrinsic 
energy splitting 2t/a, where a is the broadening of a 
single VS. We define the apparent peak separation Ae as 
the separation between the centers of the two Gaussian 
lineshapes with broadening (FWHM) a fitted to the 
calculated LDOS. The relative enhancement Ae/2t of 
the apparent splitting is most important for small t. 
For large t, the splitting of divacancy energy levels 
approaches the intrinsic energy splitting of 2t. As for 
the single vacancy level, the broadening of a divacancy 
level then scales with the square root of its relaxation 
energy (see SI). 

The hopping terms t can be extracted from Fig. 
by taking Ae equal to the observed splitting and using 
the single vacancy value 0.27 eV for a. For example, 
t = 0.11 eV was obtained from Ae = 0.28 eV for the 3NN 
divacancy. The values obtained for the hopping terms t 
are very reasonable as supported by density-functional 
theory (DFT) calculations of the unoccupied Kohn- 
Sham states of Cl divacancies in NaCl(2ML)/Cu(lll). 
As shown in Fig. |^, the calculated values for 2t of the 
various divacancy configurations are close to the values 
for the intrinsic splittings extracted from experimental 
data. However, the absolute divacancy energies are 
underestimated by the DFT calculations, as expected 
(see Fig. S3e). 

The coupling between the LISs of Cl divacancies 
and chains were studied in a simple, multiple cr-wave 
scattering model of the interface-state scattering from 
the positively charged vacancies (see methods and SI). 
The results from this model, as indicated in Fig. [^-o, 
are able to fully capture the observed behavior of the 
coupled LISs. The absence of a LIS doublet below the 


onset of the interface-state band is simply understood 
by the relatively strong interaction between the weakly 
bound LISs of the two vacancies at small intervacancy 
distances. With increasing distances between the vacan¬ 
cies, this interaction decreases, and a second LIS appears 
below the band onset. The observed symmetric charac¬ 
ter of the LIS below the band edge and antisymmetric 
character of the extended LIS above es for the 3NN 
vacancy are also revealed in the calculated LDOS images 
(see Fig. S5). Furthermore, the obtained energies for the 
LISs of the vacancy chains are in excellent agreement 
with the experiments, as indicated in Fig. 

In summary, from an AFM and STS study we find 
that the localized VSs and LISs at intentionally created 
Cl divacancies in a NaCl bilayer on Cu(lll) form sym¬ 
metric and antisymmetric states, in analogy to the bond¬ 
ing and antibonding orbitals of a hydrogen molecule. As 
expected, the energy splitting between these states in¬ 
creases with decreasing intervacancy distance. A com¬ 
parison with theory shows that the energy splitting of 
the VSs is significantly enlarged by the strong coupling 
of the tunneling electrons with phonons in the NaCl film. 
Furthermore, VSs and LISs of vacancy chains form one¬ 
dimensional quantum-well states. The VS levels of the 
chains could be well described by a simple tight-binding 
model that takes the electrostatic interaction between 
neighbouring vacancies into account. The model pro¬ 
vides also a microscopic understanding why the level 
broadening is reduced with increasing chain size. Already 
for about five coupled vacancies, a one-dimensional elec¬ 
tronic band structure evolves. 

METHODS 

Sample preparation and Cl vacancy creation. 

The experiments were carried out in a home-built low- 
temperature combined STM and AFM operated at 5K. 
A Cu(lll) single-crystal sample was cleaned by several 
sputtering and annealing cycles. NaCl was evaporated 
thermally, keeping the sample temperature at about 
270 K, such that defect-fr ee, (l OO)-terminated NaCl bi¬ 
layer islands were formecP^^. Bias voltages refer to 
the sample voltage with respect to the tip. In the 
spectroscopic measurements, the tunneling conductance 
dl /dV was recorded with conventional lock-in techniques 
with an ac bias amplitude of 25 mV at a frequency of 
294 Hz. In the double-barrier tunneling junction geome¬ 
try, the voltage drop across the insulating film will cause 
a voltage-dependent shift of the electronic levels. This 
tip-induced Stark shift is a few percent of the applied bias 
but depends only weakly on the tip distance. Therefore 
it is not considered her^Sni 

The vacancies were created 
by bringing the Cu tip into controlled contact with the 
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NaCl surface. Thereby a Cl atom is transferee! to the tip 
apex, as evidenced by a characteristic contrast change in 
STM images and a remaining depression at the prede- 
finedCl site. Constant-height AFM images with a CO 
tipj33|34| recorded to identify the vacancy location 
and exclude the generation of other close-by defects. In 
addition, we measured the local contact potential differ¬ 
ence on the polar film with Kelvin probe force microscopy 
to ensure that the defects are indeed Cl vacancieJ^. 


Theoretical calculations. 

Density functional calculations. The electronic and 
geometric structure of the divacancies and vacancy chains 
in a NaCl bilayer on a Cu(lll) surface were obtained 
using periodic density functional theory (DFT) calcula¬ 
tions. We u sed the projector augmented wave (PAW) 
methocP^I^ as implemented in VASI^. The exchange- 
correlation en ergy was described by the optB86b+vdW 
approximatiorP^*^. The incommensurate growth of 
NaCl on Cu(lll) was modeled by a super cell consist¬ 
ing of 90 Cu atoms and 32 Na and Cl atoms in each layer 
(4 layers of Cu and 2 layers of NaCl), which corresponds 
to about 2% mismatch of NaCl distances compared to 
those in bulk NaCl. 

Phonon broadening model. This simple model de¬ 
tailed in the SI contains two key parameters: the phonon 
broadening a (FWHM) of the single VS and the magni¬ 
tude of the direct interaction energy —t between the two 
VSs. The value for a is determined from the observed 
broadening of 0.27 eV of the single vacancy VS. The val¬ 
ues for t were adjusted to reproduce the observed energy 
splittings of the divacancy VSs and compared with results 
from the DFT calculations. For further information see 
the SM 

Interface-state band scattering model. The coupling 
between the LISs of Cl divacancies are modeled with a 
multiple cr-wave scattering model of the interface-state 
scattering from the vacancies. This model is detailed 
in the SP^. The interface-state band is modeled by a 
2D free-electron band with an effective mass m* and 
a band onset Eb- The values for m* = 0.46 me and 
eb = —230 meV of the interface-state band were taken 
from STM and STS measurement^^. The remaining pa¬ 
rameter eq — Eb, determines the t-matrix element for cr- 
wave scattering and was fitted to the observed value of 
-19 meV for the LIS of the vacancy. 
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FIG. 1. Vacancy states and localized interface-states of Cl divacancies, a NaCl bilayer model with a Cl vacancy 
(dashed circle) and its nearest Cl neighbours indicated, b-e Constant-height CO tip AFM images of a single Cl vacancy 
(b) and divacancies in the 5NN (c), 4NN (d) and 3NN (e) configurations, f Background-subtracted dl/dV spectrum of the 
monovacancy state (VS) peak with a Gaussian fit in red. g-j d//dV(x, V) maps along the red lines in b-e (same V scale as 
in f). Peaks in d//dV are displayed as bright, whereas the darkest areas mark regions of negative differential conductance 
(NDC)p^. The dashed line indicates the monovacancy VS peak position, k Background-subtracted d//dV spectrum of the 
localized interface-state (LIS) of a monovacancy with a Lorentzian fit in red. l-o d//dV (x, V) maps along the red lines indicated 
in b-e, respectively (same V scale as in k). The dashed line indicates the monovacancy LIS peak position and the solid line 
the interface-state band onset. White crosses indicate the vacancy positions. White arrows mark the peak positions of the 
observed d//dV resonances. Black arrows in m-o mark the peak positions of the calculated LDOS using the u-wave multiple 
scattering model. All STS measurements were acquired with Cu tips. Scale bars: 10 A. 






































FIG. 2. Vacancy states and localized interface-states 
of Cl vacancy chains, a Constant-height CO tip AFM 
images of Cl vacancy chains with N = 4, 5, 6 vacancies in 
the 5NN configuration (apolar direction), b Corresponding 
d//dV(x,V) maps of the LISs along the chains. The con¬ 
tinuous line indicates the interface-state band onset. White 
arrows mark the peak positions of the observed dl/dV res¬ 
onances, and black arrows mark the peak positions of the 
calculated LDOS. c Constant-height CO tip AFM images of 
Cl vacancy chains with V = 3, 4, 5 vacancies in the 3NN con¬ 
figuration (polar direction), d Corresponding dl/dV{x,V) 
maps of the VSs along the chains. White crosses mark the 
positions of the vacancies and the white arrows dl/dV peaks. 
Scale bars: 10 A. 
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FIG. 3. Tight-binding models, a Experimental VS broad¬ 
ening fj as a function of the 3NN chain length, b Tight- 
binding (TB) model for a V = 3 chain, c TB model in¬ 
cluding the electrostatic on-site energy shift between nearest 
neighbours, d Experimental VS energies as a function of the 
3NN chain length, e TB model energies obtained by fitting 
the on-site energy Eq = 2.8 eV and hopping term t = 0.14eV 
to the experimental mono- and divacancy VS energies shown 
in d. f TB model accounting for a shift of the on-site energy 
{S = 0.07 eV) due to the electrostatic interaction between 
nearest neighbours. 
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FIG. 4. Lineshapes of vacancy states of a Cl diva¬ 
cancy. a,b Calculated adiabatic potential energy surfaces 
(PESs) of the electronic ground-state (black) and the elec¬ 
tronically excited and singly occupied divacancy states (red 
and blue) as a function of symmetric qs (a) and antisymmetric 
Qas (b) phonon displacements. The vibrational ground-state 
probability distribution functions (pdf) of Qs (a) and Qas (b) 
are indicated by the green curves. Here, 2t is the intrin¬ 
sic electronic level splitting. The inset in (a) shows Qs and 
Qas in terms of the displacements qi and q 2 of the two Ein¬ 
stein modes each associated with a vacancy. The black crosses 
mark eqnilibrium displacements of vacancy 1 and 2 when the 
divacancy is occupied by a single electron. c,d Energy-level 
diagram for the divacancy states as a function of symmetric 
Sss (c) and antisymmetric Ssas fluctuations (d) of the mono¬ 
vacancy energies induced by the phonon motion. The vertical 
gray arrows indicate the electron-attachment process and the 
black curves the resulting LDOS. e Calculated LDOS of the 
3NN divacancy for uncorrelated vacancy energy fluctnations 
(red) compared to a superposition of two Gaussians separated 
by 2t and each with a broadening a (blue). Ae denotes the 
level splitting for nncorrelated vacancy energy fluctnations. f 
Scaled level splitting Ae/cr as obtained from a Gaussian fit 
to the calculated lineshapes of the divacancy states for un¬ 
correlated vacancy energy fluctuations (red) compared to the 
intrinsic splitting (bine) as a fnnction of the scaled intrinsic 
level splitting 2t/a. The calculated and observed scaled split¬ 
tings for different divacancy configurations are indicated by 
black squares and dashed horizontal lines, respectively. The 
green curve shows the apparent level broadening a for uncor¬ 
related vacancy energy fluctuations. 













































